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Comment on "Lyapunov Exponent of a Many 
Body System and Its Transport Coefficients" 

In a recent Letter, Barnett, Tajima, Nishihara, Ucshima 
and Furukawa [EI obtained a theoretical expression for 
the maximum Lyapunov exponent Ai of a dilute gas. 
They conclude that Ai is proportional to the cube root of 
the self-diffusion coefficient D, independent of the range 
of the interaction potential. They validate their con- 
jecture with numerical data for a dense one-component 
plasma, a system with long-range forces. We claim that 
their result is highly non-generic. We show in the fol- 
lowing that it does not apply to a gas of hard spheres, 
neither in the dilute nor in the dense phase. 

Systems of hard spheres have properties similar to real 
fluids and solids and provide a reference for successful 
perturbation theories [Q. Simulations with this model 
were able to uncover fundamental aspects of collective 
particle dynamics such as recollisions and the "cage" 
effect Hard-sphere systems are also paradigms for 
the chaotic and ergodic properties of many-body systems 
with short range interactions, and were shown to have a 
positive Kolmogorov-Sinai entropy ||[|. 

For dilute gases, Krylov H provided an analytical esti- 
mate for the maximum Lyapunov exponent, 


(32nK/3mN) 1/2 a 2 n\og (na 3 n/V2 


(1) 


where K is the kinetic energy, N is the number of par- 
ticles, to the particle mass, n the number denstiy, and 
a is the hard sphere diameter. This expression has been 
verified numerically (apart from a factor ~ 2.8 Q), and 
has been extended to larger densities 0] . 

The diffusion coefficient for dilute hard-sphere gases is 
well approximated by the Enskog expression 
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A comparison of Eqs. ([l]) and (g) reveals that, in the 
dilute gas limit, the proposed relation Ai cx D 1 / 3 of Bar- 
nett et al. cannot be satisfied. Moreover, we combine in 
Fig. 1 recent simulation results for D and Ai, which were 
obtained for a system of 500 hard spheres over the full 
range of fluid densities (0.0001 < na 3 < 0.89). Reduced 
units are used for which er, to, and the kinetic energy per 
particle K/N are all unity. One observes that these data 
are not consistent with the proposed _D 1 / 3 -dependence 
(solid line), neither for low densities nor for large. 

We conclude that the conjecture by Barnett et al. does 
not apply to many-body systems with short-range inter- 
actions. But even its applicability for long-range inter- 
actions is doubtful. A one-dimensional gravitational sys- 
tem with finite N exhibits a positive Ai [g, whereas this 
clustering and confining system does not show diffusion. 
We also note that, while the theoretical expression (26) 


in Ref. M has been obtained for a dilute gas, the data in 
Fig. 1 of Ref. § are for a dense plasma with a Coulomb 
coupling costant T ranging from 1 to 150. As reported by 
the same authors || , for T > 1 the plasma behaves as a 
liquid and not as a gas. The dilute gas limit is recovered 
only for T< I §. 

We thank M. Antoni, U. Balucani, A. Rapisarda, and S. 
Ruffo for useful discussions. 
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Fig.l: Simulation results Ai = Ai(Z?) (crosses) for a gas of 
hard spheres. The solid line refers to the expression Ai CX D 1 / 3 
suggested in Ref. Jl| 
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